// 利用排列组合计算所有的可能情况的数量 365^n
// 计算生日不在同一天的情况的数量 C_365^n
// 最终的可能性为 1 - C_365^n/365^n
// 化简为 1 - 365-1/365 * 365-2/365 * ...... * 365-n+1/365
pub fn new_birthday_probability(n: u32) -> f64 {
    let mut p = 1.;
    for i in 1..n {
        p *= (365_f64 - i as f64) / 365_f64;
    }
    1. - p
}
